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ABSTRACT: We introduce a new measurable quantity, i?A</>> f° r studies of the rapidity and 
transverse momentum dependence of dijet azimuthal decorrelations in hadron-hadron colli- 
sions. In pQCD, i?A</> is computed as a ratio of three-jet and dijet cross sections in which the 
parton distribution functions cancel to a large extent. At the leading order, is propor- 
tional to a s , and the transverse momentum dependence of i?A0 can therefore be exploited 
to determine a s . We compute the NLO pQCD theory predictions and non-perturbative 
corrections for i?A</> at the LHC and the Tevatron and investigate the corresponding uncer- 
tainties. From this, we estimate the theory uncertainties for a s determinations based on 
Ra</> at both colliders. The potential of i?A0 measurements for tuning Monte Carlo event 
generators is also demonstrated. 
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1 Introduction 

Theory predictions for inclusive jet and dijet cross sections in hadron-hadron collisions at 
fixed order in perturbative Quantum Chromodynamics (pQCD) are currently available at 
next-to- leading order (NLO) in the strong coupling constant a s . Using precise experimental 
data, these predictions have been well tested and applied in determining the parton distri- 
bution functions (PDFs) of the proton, and a s [1]. Direct tests of pQCD at higher orders 
require measurements of quantities probing multi-jet final states with three or more jets. 
Quantities in which a cross section for the production of three or more jets is normalized 
by a dijet cross section (or an inclusive jet cross section) are ideal for a s determinations. 
These quantities are still sensitive to the degrees of freedom in the multi-jet final state and, 
in pQCD, proportional to (at least) 0(a s ), while the PDF sensitivity exhibited by a typical 
multi-jet cross section [2] can be strongly reduced. Examples of such quantities are the ratio 
of the inclusive three jet and dijet cross sections, R3/2 [3-5], and the average number of 
neighboring jets, Rar, which has recently been proposed, measured, and used to determine 
a s [6]. 

A third related quantity is the dijet azimuthal decorrelation, which studies the relative 
angle in the azimuthal plane between the two jets with the highest transverse momentum 
(pr) A<^dijct = 1 0jeti — 0jet2 1 • hi calculations at 0(a 2 s ) , dijet events have exactly two jets with 
equal pt, and their azimuthal angles are correlated such that A^ijet = ft- Deviations from 
A^dijct = 7r (hereafter referred to as "azimuthal decorrelations") are caused by additional 
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Figure 1. A sketch of the angle A^dijet in the azimuthal plane, in dijet events for different amounts 
of additional radiation outside the dijet system. 

radiation which is not clustered into the two highest px jets. Additional radiation with 
small pt reduces A(fidijet by a small amount, while high-px radiation can reduce A^ijet 
significantly thereby leading to larger azimuthal decorrelations as illustrated in figure 1. 
Due to kinematic constraints, three-jet final states are restricted to A^ijct > 2-7r/3, while 
the phase space of Ac^dijet < 2-7r/3 is only accessible in final states with at least four jets. 

The D0 collaboration has introduced the quantity (l/a"dijet) ■ ^Odijet/rfA^dijeti which 
is the dijet cross section differentially in A^ijetj normalized by the inclusive dijet cross 
section Odijet ( m the same kinematic range and integrated over A0dij e t) [7]. This quantity 
was measured in pp collisions at y/s = 1.96 TeV, for different p™ ax requirements, where p™ ax 
is the highest jet px in the event, and for a fixed px requirement for the second leading px 
jet. For this quantity, the range from small to large azimuthal decorrelations can be used to 
study the transition from soft to hard higher-order pQCD processes and the measurement 
results placed strong constraints on Monte Carlo parameters [8]. The same analysis strategy 
was later employed by the CMS and ATLAS collaborations using pp collision data at -y/i = 
7 TeV, thus accessing larger p™ ax [9, 10]. The common approach focuses on the A^ijet 
dependence; the px dependence is not easily visible in these presentations. Furthermore, in 
pQCD, dijet azimuthal decorrelations are predicted to depend not only on px, but also on 
the rapidities of the two leading px jets. The measurements by the D0, CMS, and ATLAS 
collaborations, however, did not explore the rapidity dependence. 

In this article, we propose a new quantity for studying dijet azimuthal decorrela- 
tions with emphasis on the rapidity and the px dependence. 1 The former aspect will allow 
us to perform novel tests of the pQCD predictions, while the latter can be exploited for 
determinations of a s and its running. The article is structured as follows: In section 2 
we motivate the variables used to study the rapidity and px dependencies, and we give 
the definition of R^. In addition, we propose realistic scenarios for phase space regions in 
which i?A^ can be measured by the LHC and the Tevatron experiments. Theory predictions 
for these scenarios are presented in section 3, including perturbative and non-perturbative 
contributions. The possible impact of R^ measurements on determinations of a s and on 
Monte Carlo tuning is discussed in section 4. 

1 Some initial studies for experimental measurements of i?A0 have been made in references [11, 12]. 
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Figure 2. Jet rapidity variables in the transverse plane in the dijet center-of-mass frame and in 
the hadron-hadron center-of-mass frame. 

2 Definition and Phase Space Scenarios 

The quantities -R3/2 [3-5] and Rar [6] are defined as ratios of multi-jet cross sections. These 
ratios can be interpreted as the conditional probability that an event with two high-p^ jets 
also contains a third jet (-R3/2) an d as the average number of neighboring jets for a given 
jet (Rar)- We propose to study dijet azimuthal decorrelations, using a quantity with a 
similar intuitive interpretation. For this purpose we introduce the quantity Ras- Before 
we define Ra<p, we motivate the variables used to study the rapidity and pt dependencies. 
In addition, we propose realistic scenarios for measurements of Ra</> at the LHC and the 
Tevatron. 

2.1 Variables for the Rapidity and pt Dependence 

The Rapidity Variable One of the main goals for the new quantity Ra4> is to measure 
the rapidity dependence of dijet azimuthal decorrelations. The previous analyses [7, 9, 10] 
applied rapidity requirements for both jets in the hadron-hadron center-of-mass frame (i.e. 
the lab frame). In general, this frame is, however, longitudinally boosted with respect 
to the center-of-mass frame of the hard subprocess (corresponding to yboost) ; as shown in 
figure 2. In the approximation of 2 — > 2 processes, the rapidities y\ and y\ (in the dijet 
center-of-mass frame) have the same magnitude (y* = |y*| = {y^D, and are related to the 
rapidities y± and 2/2 (in the hadron-hadron center-of-mass frame) by 

yi=y* + 2/boost and y 2 = y* + y hoost . (2.1) 

We propose to measure the rapidity dependence of Ra4> as a function of the variable y* for 
a fixed requirement for the variable yboost • 

Both variables are given by 
2/boost = (yi + 2/2)/2 and y* = \y x - y 2 \/2 , (2.2) 
where y\ and y 2 are the respective rapidities of the two leading pj- jets in the event. 

The Pt Variable In the leading logarithmic approximation, an n-parton final state can 
be regarded as emerging from a two-parton final state through successive branching, as 
displayed in figure 3 for three- and four-parton final states. In this picture, a quantity such 
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Pn = Pit Pn < Pn' 

(Pn + Pt 2 ) < (Pit + Pt?) (Pn + Pt 2 ) < (Ptv + Pre-) 

H T = Ej (p Ti ) - E r (p Ti .) = H' T H T = E, (p Ti ) - E r (p Tr ) = H' T 



Figure 3. Sketches of the azimuthal plane in which a three-parton (a) and a four-parton final state 
(b) emerge from a two-parton final state through parton-branching, and the relations of different 
variables before and after the branching. The variable Ht is approximately preserved in these 
processes. 

as -R3/2 may be interpreted as the branching probability from two to three final state partons 
(for partons which have a sufficient angular separation to be resolved as individual jets, 
according to the jet definition). However, this interpretation only holds if the quantity is 
binned in a "p^-type" variable that does not change its value before and after the branching 
of the third parton. Examples of such variables are the leading jet pr, or, approximately, 
Ht (defined as the scalar pt sum over all jets in an event). A counter example is the 

(2) 

variable H T = (pri + PT2), defined as the scalar pt sum of the two leading jets, which 
is reduced after the branching displayed in figure 3 (a). For a quantity like the ratio of 
inclusive four-jet and dijet cross sections, R4/2, the leading jet pt could also be reduced, 
e.g. by a branching as displayed in figure 3 (b). 

In the measurement of dijet azimuthal decorrelations, the angle A0dij e t is affected not 

only by the third jet but by all additional radiation in the event. Therefore variables like 

(2) 

p™ ax or H T are not suited if we want to interpret the ratio i?A</> as the probability for parton 
branching. Only the value of Ht is approximately conserved after the branching processes. 
With an ideal detector and in a clean environment, one might want to define Ht as the 
scalar pt sum of all jets in the event, without any pt or y requirements. In practice, one 
has a limited detector y acceptance, a limited knowledge of the detector response for low pt 
jets, plus contributions from the underlying event. Therefore it is advisable to restrict the 
Pt sum in the Ht definition to jets which are well measured, and for which non-perturbative 
contributions are small, by requiring pTi > PTmm- A limited detector y acceptance can be 
taken into account by requiring that the jets are contained inside this acceptance region. 
Since we study the rapidity dependence based on y* which is longitudinally boost invariant, 
we would like to preserve this property also for the Ht definition. Therefore the jet selection 
is not based on the absolute jet rapidities in the lab frame, but on the longitudinally 
boost invariant quantity — yboostl, an d Ht is computed as 

H T = Y,PTi, (2.3) 

iec 
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based on all jets in the set C which is defined as 

C = { i | 1 < i < rijet; and p Ti > p T mm, and \y { - y hoost \ < y max } . (2.4) 

In this definition, nj e t is the total number of jets in the event, and prmin an d y max are 
parameters which can be chosen according to the experimental environment. The value of 
y max should be chosen at least as large as the maximum accessible y* for two leading pt 
jets, to ensure that these are always members of the set C , and therefore included in the 
Ht sum. 

2.2 Definition of 

With the criteria above, we propose to study dijet azimuthal decorrelations using the new 
quantity i?A$, which represents the fraction of all inclusive dijet events for which the two 
leading px jets have a decorrelation of A</>dijet < A^ max . It is defined as 

R A ,(H T ,y*,A^ x ) = dVd ^; sive) . (2.5) 

dHj- dy* 

The denominator, cPudijet (inclusive) / (cIHt dy*), is the cross section for the production of 
two or more jets, with pt > PTmm, an d 2/boost < 2/boosf double differentially in the variables 
y* and Ht- The numerator, <i 2 <Tdij C t(A</>dijct < A^> max )/(<iffT dy*), is a subset of the denom- 
inator with the additional requirement that the two leading pt jets have A^dijct < A0 max . 
The quantity is measured as a function of the parameter A</> max , and in bins of y* and 
Ht, and therefore expressed as Ra4>(Ht,U*, A</> max ). 

It may be convenient to introduce an additional requirement of an if^-dependent lower 
limit on the leading jet pt as pri > f ■ Ht, in both the numerator and the denominator. 
This requirement (which cuts the tail of low leading jet pt) is necessary in the experiment 
if events are triggered by inclusive single jet triggers. The value of / should not be too 
large, so as to not restrict the multi-jet phase space too strongly. We recommend to set 
/ not larger than / = 1/3, so that the phase space for 2 — > 2 and 2 — > 3 processes is not 
affected. 

2.3 Phase Space Scenarios for the LHC and the Tevatron 

To produce specific theory predictions, we propose two scenarios of phase space regions in 
which Ra<i> can be measured at the LHC and the Tevatron. While making realistic choices 
that take into account current practices by the experiments, we try to keep the two scenarios 
as similar as possible, so that the results can be used to study the y/s dependence of i?A^- 

LHC Scenario We assume the running conditions of 2012, where the LHC was producing 
pp collisions at y/s = 8TeV. Following the choices by the ATLAS and CMS experiments, 
jets are defined using the anti-kj jet algorithm [13], here with a jet radius of R = 0.6 (in 
the y-(p plane), which is within the range of 0.4-0.7 that is used by ATLAS and CMS. The 
parameters in the R^ definition are set to y max = 2.0, y™^ st = 0.5, and PTmm = 100 GeV. 
The additional iJ^-dependent requirement on the leading jet pt is pti/Ht > 1/3. The 
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nun 11 cone, K conc =0.7 


PTmin (.^TminJ 


100 GeV (0.0250) 


30 GeV 0.0306) 


max 
"boost 


O.o 


O.o 


* 

Umax 


2.0 


2.0 


Pti/Ht 


> 1/3 


> 1/3 




7vr/8, 5vr/6, 3vr/4 


7vr/8, 5vr/6, 3vr/4 


y* ranges 


0.0-0.5, 0.5-1.0, 1.0-2.0 


0.0-0.5, 0.5-1.0, 1.0-2.0 


Ht range 


750-4000 GeV 


180-900 GeV 



Table 1. Summary of the phase space definitions for the LHC and the Tevatron scenarios. 

Umax an d yboost requirements ensure that the two leading pt jets and all other jets entering 
the Ht sum are well-contained in the detector, within \y\ < 2.5. The PTmin requirement 
ensures that all jets are well measured in the experiment, and that pileup contributions and 
non-perturbative corrections are small. For studies of the -y/s dependence of R/\(f>, we note 
that the pTmin requirement translates to a requirement for the scaling variable xt = Ip^j \fs 
of XTmin = 0.025. The parameter A0 max is set to 77r/8, 5-7t/6, or 37r/4, and the y* regions 
are chosen as0<y* <0.5, 0.5<y* < 1. and 1 < y* < 2. The Ht dependence is studied 
over the range 750 < H T < 4000 GeV. 

Tevatron Scenario In Run II, the Tevatron collided protons and anti-protons at y/s = 
1.96 TeV. For the majority of the jet results, the CDF and D0 experiments use iterative 
seed-based cone algorithms with a cone of radius R CO nc = 0.7 in y and (j). For these studies, 
we apply the Run II midpoint cone jet algorithm [14] that is used by D0. We use the same 
values for the parameters y max , y^oo St , an d A^ max , the identical y* regions, and the same 
Pti/Ht requirement as in the LHC scenario. The only differences are the value of the PTmin 
requirement, which is set to prmin = 30 GeV, and the Ht range of 180-900 GeV. The prmin 
requirement translates to a requirement for the scaling variable of XTmin ~ 0.0306, which is 
slightly higher than the corresponding requirement in the LHC scenario (of XTmin = 0.025). 
However, we use this value because it corresponds to the lower pt requirements used in 
recent multi-jet measurements at the Tevatron [5, 6] 

The parameters, defining the phase space for the LHC and the Tevatron scenarios are 
summarized in table 1. 

3 Theory Predictions and their Properties 

In this section, we compute the perturbative and the non-perturbative contributions for 
i?A0 and investigate their properties. We compare the predictions for the LHC and the 
Tevatron scenarios, and investigate the differences due to the differences in -y/i, and due to 
the slightly different phase space requirements. 

All theory results have been obtained using the implementations of the anti-k^ and the 
D0 Run II cone jet algorithms in FASTJET [15. 16]. 
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Figure 4. The NLO pQCD predictions for Rac/> as a function of (Ht/^/s), and in regions of 
y* (columns) and for different A0 max requirements. The results are shown for the LHC and the 
Tevatron scenarios, and for a LHC-like scenario in which the LHC scenario is scaled to ^/s = 2 TeV 
and modified to pp collisions. 



3.1 NLO pQCD Predictions 

The NLO (LO) pQCD predictions for Ra6 are computed by taking the ratios of the 
NLO (LO) pQCD predictions for the cross sections in the numerator and the denomi- 
nator in eq. (2.5). The denominator is the inclusive dijet cross section for which the NLO 
(LO) prediction is computed at O(a^) (0(a 2 )). Due to the additional requirement of 
A(/>dijet < A</> max , the numerator receives only contributions from final states with three or 
more jets. Therefore, the numerator is a three-jet cross section for which the NLO (LO) 
prediction is computed at O(ctf) (0(a^)). All NLO and LO pQCD results are computed 
using NLOJET++ [17, 18], interfaced to FASTNLO [19]. The calculations are made in the MS 
scheme [20] for five active quark flavors, and using the next-to- leading logarithmic (two- 
loop) approximation of the renormalization group equation. The value of a s (Mz) = 0.118 
is used consistently in the matrix elements and in the MSTW2008NLO PDF sets [21]. The 
central choice /j,q for the renormalization and factorization scales is [ir = fip = fiQ = Ht/2. 2 
The results of the NLO calculations for the LHC and the Tevatron scenarios are dis- 
played in figure 4, where R&a, is shown as a function of (Ht/\J~s) in different regions of y* 
and for different A^ max requirements. In different regions of Hf and y* , and for different 
choices of A(/> max , has values in the range 0.012-0.32. In most phase space regions, 

i?A0 decreases with increasing Ht, except at 1 < y* < 2 where i?A</> increases again at high 

2 At LO, this choice coincides with the common choices of [ir.f = Pt for inclusive jet production and 
fiR,F = (pti +PT2)/2 for dijet production. 
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Ht- At fixed Ht, Ra^ increases with increasing y*. The fact that Rac/> decreases with 
decreasing A0 max is a trivial phase space effect, since a stronger Ac^ijet requirement leads 
to a smaller cross section in the numerator. 

For a fixed XTmin requirement and at fixed (Ht/>/s), the y/s dependence of the per- 
turbative results for Rai/> is only introduced through the evolution of a s and the PDFs with 
the scales /i/j and //p. In figure 4, the y/s dependence of Ra4> cannot directly be judged 
based on the comparison of the LHC and the Tevatron scenarios, as the two differ in the 
a?Tmin requirement and in the jet algorithm. The following study is made to separate the 
latter effects from the genuine y/s dependence of Ra6- Using the flexibility provided by 
FASTNLO, we use the FASTNLO coefficient tables for the LHC scenario (for pp collisions at 
y/s = 8TeV) to compute the corresponding predictions for pp collisions at the same y/s. 
The results for the latter (not shown in figure 4) agree with those for the LHC scenario 
better than 0.8% for Ht < 2TeV and always better than 3.2% in the phase space studied, 
meaning that Ra^> is insensitive to the difference between pp and pp initial states. Then we 
use the FASTNLO results for the LHC scenario to compute the corresponding predictions for 
a LHC-like scenario (i.e. using the same jet algorithm and the same XTmin requirement) for 
pp collisions at y/s = 2 TeV. These predictions are shown in figure 4 as the dotted line. The 
results at y/s = 2 TeV are 10-20% higher than those at 8 TeV. This y/s dependence is 
consistent with the running of a s over a factor of four in energy. 3 For A<^ max = 7tt/8 and 
5tt/6, the results for the LHC-like scenario agree within 5% with those for the Tevatron 
scenario. Only for A0 max = 37r/4 the differences become larger (slightly more than 10% 
at y* > 1). From this we conclude that, even with different jet algorithms and slightly 
different XTmin requirements, a comparison of Tevatron and LHC data is probing the y/s 
dependence of Ra<)> and testing the corresponding theory predictions. 

In the following, we investigate the NLO /c-factors and the scale dependence as indica- 
tors for the stability of the perturbative expansion, and we study the PDF uncertainties for 
Ra6- The NLO /c-factors are computed as the ratio of the NLO and the LO predictions, 
k = Ra^/RaI- The values of the & - f act 

ors are displayed in figure 5 as a function of 
(Hx/y/s), for the LHC and the Tevatron scenarios. For A0 max = 7ir/8 and 5n/6, the k- 
factors for the LHC and the Tevatron are always close to unity; they decrease slightly with 
increasing y* and are almost independent of Ht- Due to kinematic constraints, the region 
of A^dijct < 27r/3 is only accessible in four-jet final states. For this reason, the kinematic 
region of A^ max = 37r/4 also receives large contributions from four-jet production which are 
only modeled at LO by the 0(a^ ) calculation for the numerator of Raa- This is reflected 
in the large NLO k- factors for A</> max = 37r/4 which are as large as k = 1.5 at lower Ht- 

The uncertainties due to the scale dependence are computed from the relative variations 
of the RA(f> results when /i/j and \ip are varied independently around /j,q = Ht/2 between 
fj,o/2 and 2//o but never exceeding 0.5 < Hr/hf < 2.0. These uncertainties are displayed 
in figure 6 for the LHC and Tevatron scenarios. For the LHC (Tevatron) scenario, these 
uncertainties are typically 5% (7%) for A<fi max = 7ir/8 and 57r/6, and up to 19% (17%) 

3 The PDFs approximately cancel in the ratio Ra4>, so the hf dependence of the PDFs does not have a 
significant impact on the y/s dependence of Ra<i>- 
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Figure 5. The NLO fc-factors for Ra<p as function of (Ht/^/s), in different regions of y* (columns) 
and for different A</> max (rows), for the LHC and the Tevatron scenarios. 



for A(p 

max — 37r/4. The latter is directly related to the large NLO fc-factors in this region. 
The uncertainties for A(/> max = 7n/8 and 57r/6 are slightly smaller as compared to the 
uncertainties for other ratios of three-jet and dijet cross sections like R3/2 and Rar as 
recently measured at the LHC [4] and the Tevatron [5. 6]. 

The PDF uncertainties are computed for the MSTW2008NLO PDFs using the up and 
down variations of the 20 orthogonal PDF uncertainty eigenvectors, corresponding to the 
68% C.L. The MSTW2008NLO PDF uncertainties for the LHC and Tevatron scenarios are 
shown in figure 6 and they are always below 1%. Also displayed in figure 6 are the NLO 
pQCD predictions obtained for CT10 [22] and NNPDFv2.1 [23] PDFs. Those results are 
larger by up to 5% as compared to the results for MSTW2008NLO PDFs. The largest 
deviations occur at smallest A^ max . 

From these studies we conclude that theory predictions are most reliable (as indicated 
by a small scale dependence and £>factors which are close to unity) in the kinematic regions 
of y* < 1 and for A<^> max = 7tt/8 and 57r/6. 

3.2 Non-perturbative Effects 

In these studies, we consider non-perturbative effects due to the underlying event and due to 
hadronization corrections. Both corrections are estimated using the models implemented in 
the event generators HERWIG 6.520 [24, 25] and PYTHIA 6.426 [26, 27]. The HERWIG results 
are obtained using default settings and the PYTHIA results are obtained for four different 
popular tunes. These are tunes DW [28] and A [29], which use a Q 2 ordered parton shower 
and an older model for the underlying event, and the tunes AMBT1 [30] and S Global [31], 
which use a pt ordered parton shower and the new model for the underlying event [32, 33]. 
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Figure 6. The renormalization and factorization scale dependence and the PDF uncertainties for 
the MSTW2008NLO PDFs for the NLO pQCD predictions for the LHC and Tevatron scenarios. 



The non-perturbative corrections are obtained from the results of three calculations in which 
R&4> is computed 

(1) at the parton level (the partons after the parton shower) with no underlying event, 

(2) at the particle level (using all stable particles) with no underlying event, and 

(3) at the particle level (using all stable particles) with underlying event. 

The total non-perturbative correction c npcr t is defined as the product of the hadronization 
correction Chadr an d the underlying event correction c uo which are each given by ratios of 
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Figure 7. The hadronization corrections for Ra4>, plotted as a function of Ht in different y* 
regions (columns) and for different values of AcA max (rows), for the LHC and Tevatron scenarios. 



Ras results on different levels, as 



R 



Cnpert — Cli a clr ' Cv 



(2) 
Acf> 



R 



(3) 
A<f> 



with c hadr = — j^- and c ue = — . (3.1) 

Figure 7 shows the hadronization corrections for Ra4> for the LHC and the Tevatron 
scenarios. The hadronization corrections for the LHC are very small and always be- 
low 1.5% (0.985 < Chadr < 1-00) at all Ht-, y* , and for all At/> max requirements. The 
hadronization corrections for the Tevatron, although slightly larger, are still always below 
6% (0.94 < Chadr < 1-01). The HERWIG results and the PYTHIA results for the different 
tunes agree always within 1% (3%), for the LHC (Tevatron) scenario. 
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Figure 8. The underlying event corrections for Rac/>, plotted as a function of Ht in different y* 
regions (columns) and for different values of Ac/> max (rows), for the LHC and Tevatron scenarios. 



The corrections for i?Ac6 due to effects from the underlying event are displayed in figure 8 
for the LHC and the Tevatron scenarios. For the LHC, the underlying event corrections are 
always below 4% (0.99 < c uc < 1.04), at all H Tl y* , and for all Ac/> 

max requirements. The 

maximum corrections increase with decreasing A(/> max . They are 2% for A(/> max = 7n/8, 
3% for Aci max = 57r/6, and 4% for Ac/> max = 3tt/4. For the Tevatron, the underlying 
event corrections (0.99 < c uc < 1.08) have the same qualitative behavior and they are 
approximately twice as large as those for the LHC. The HERWIG results and the PYTHIA 
results for the different tunes are in good agreement. The different model predictions for 
the LHC (Tevatron) scenario agree better than 2% (3%) for Ar/> max = 7n/8, and always 
better than 4% (6%). 
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The total non-perturbative corrections are always in the range 0.98 < c npcr t < 1.03 
(0.96 < c npcrt < 1.03) for the LHC (Tevatron) scenario. The smallness of these corrections 
and their small model dependence are remarkable features of the quantity which makes 
RAcf> well-suited for precision tests of pQCD. 

4 Phenomenology 

In this section, we discuss two examples of the potential impact of future i?A</> measurements 
for QCD phenomenology. 

4.1 Sensitivity to a s 

In pQCD, Ra<p is computed as a ratio of three-jet and dijet cross sections, which is, at LO, 
proportional to a s . In the following, we study the sensitivity of R^ to a s and investigate 
the effects of theoretical and experimental uncertainties on the a s results. These studies 
are made in the kinematic region of y* < 0.5 and for A0 max = 7ir/8 where the pQCD 
predictions are most reliable (see section 3.1). 

In all studies, a s is varied consistently in the pQCD matrix elements and in the PDF 
sets. The continuous dependence of the NLO pQCD predictions for R^ on a s is obtained 
using cubic interpolation between the discrete a s (Mz) values for which the MSTW2008NLO 
PDFs sets are available. Where needed, a s (Mz) is converted from the scale [Ir = Mz to 
the scale ^jlr = Ht/2, using the two-loop solution of the renormalization group equation. 

A first impression of the a s sensitivity is obtained by studying the a s (Mz) depen- 
dence of the NLO pQCD predictions for R^4>- For this purpose, we plot i?A</> for values 
of a s (Mz) = 0.110-0.130 (labeled R&<p{a s {Mz)))-, normalized to the value of i?A0 for 
a s (Mz) = 0.1184 (the world average value [1]) as a function of a s {Mz)- The results are 
shown in figure 9 for three different Ht bins, for the LHC and the Tevatron scenarios, and 
are compared to the naive expectation of a linear relation (Ra6 (X a s ). Deviations from a 
linear dependence could be due to three different effects. 

1. The naive expectation of a linear dependence stems from the LO picture, and is 
modified due to NLO corrections. 

2. The naive expectation assumes a perfect cancellation of the PDFs, while residual PDF 
effects may lead to non-linearities. 

3. While figure 9 shows the a s (Mz) dependence, the calculations for i?A</> are made for 
the scale [1r = Ht/2, and the relation between a s (fj,ji) and a s (Mz) as a function of 
a s (Mz) is not linear, and involving logarithms of {^r/Mz)- 

For a s (Mz) < 0.125, in the Tevatron scenario, the a s (Mz) dependence of Ra</> is almost 
linear for 180 < Ht < 205 GeV (i.e. where [Ir = Ht/2 ~ Mz)- The decrease of the slope 
(and therefore the increasing non-linearity) of the curves for higher Ht is likely caused by 
the third effect. The change of the slopes around a s (Mz) ~ 0.125 is likely caused by the 

4 The MSTW2008NLO PDF sets are available for a s (Mz) = 0.110, 0.111, 0.112, • • • , 0.130. 
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Figure 9. The a s (Mz) dependence of Ra</>, normalized by the Ra<p value at a s (Mz) = 0.1184, for 
the LHC scenario (left) and the Tevatron scenario (right) in three different Ht bins for A0 max = 
77r/8 and if < 0.5. 



second effect. The reduced slope towards high Ht implies a slightly reduced sensitivity to 
a s (M z ) at the LHC. 

Currently, the precision of a s results obtained from hadron colliders which are based 
on NLO calculations, is limited by theory uncertainties stemming from the renormalization 
and factorization scale dependencies of the calculations. 5 Therefore we estimate the corre- 
sponding uncertainties for a s extractions from R/\^. In the typical procedure of most a s 
analyses, the central a s results are derived for a fixed choice of the renormalization and 
factorization scales. The uncertainties of a s due to the scale dependence are then obtained 
by repeating the a s fits for variations of the scales around their central values. In the 
absence of actual data, we estimate the corresponding uncertainties for a s by comput- 
ing the variations in a s which are required to bring the NLO pQCD results at a different 
scale into agreement with those at the central scale. As discussed in section 3.1, we use 
a central scale of [Ir = /ip = fio = Ht/2 and a range of variations in which /ir and \xp 
are varied independently between fio/2 and 2//q, while never exceeding 0.5 < ^r/hf < 2.0. 
The largest effects of all variations are quoted as the corresponding uncertainties for a s . 
The expected uncertainties for 0£ s (fJ,R = Hf/2), derived using this procedure are shown in 
figure 10 (left) in the kinematic range y* < 0.5 and A<r5 max = 7ir/8, as a function of /ir, 
for the LHC and the Tevatron. For the LHC scenario, the uncertainties are between —3% 
and +4% at high Ht, and slightly lower at low Ht- For the Tevatron, the uncertainties 
are between —4% and +6%. This is a good theoretical precision for testing the running of 
a s at highest energies. 

5 The presently most precise a B result from a hadron collider was obtained using theory calculations 
beyond NLO (adding the 2-loop corrections from threshold corrections) and has therefore smaller scale 
uncertainties [34]. These contributions are, however, only available for inclusive jet production and neither 
for dijet nor for three-jet production. 
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Figure 10. The estimated uncertainties of a s results extracted from Ra<p at the LHC and the 
Tevatron, at a scale /in = Ht/2, due to the renormalization and factorization scale dependence of 
the NLO pQCD calculation (left) and due to an experimental uncertainty of 4%, for Ac/> max = 7tt/8 
and y* < 0.5. 



For comparison, we have also computed the uncertainty of oi s {Ht/2) resulting from an 
experimental uncertainty of 4%. The results in figure 10 (right) show that for the LHC and 
the Tevatron scenarios this uncertainty is approximately of the same size as the uncertainty 
due to the scale dependence. In other words, if the total uncertainty for the a s results 
should not be limited by the experimental precision, the experiments must measure R/\c/> 
with a precision of at least approximately 4%. Given the precision of recently published 
measurements of the multi-jet cross section ratios R3/2 [3-5] and Rar [6], this should be 
achievable. 



4.2 Event Generator Tuning 

We also compute R^ predictions using the Monte Carlo event generators HERWIG 6.520 and 
PYTHIA 6.426. The HERWIG generator is used with default settings, and for PYTHIA we use 
tune AMBT1 (derived by the ATLAS collaboration using LHC data) and tune DW (which 
has been tuned to the previous D0 measurement of dijet azimuthal decorrelations) . The 
results for Ac/> max = 57r/6 are shown in figure 11 for the LHC and the Tevatron scenarios. 
At central rapidities (i.e. low y*) the predictions of the different generators agree reasonably 
well with each other, and also with the NLO pQCD results. The different generators (and 
tunes), however, predict very different y* dependencies. In the forward region (1 < y* < 2), 
the predictions differ strongly in magnitude and in shape. The HERWIG predictions are more 
than 30% higher, and the PYTHIA tune AMBT1 predictions have a very different shape as 
compared to the NLO pQCD results. The large range of the predictions from different 
PYTHIA tunes and their differences to HERWIG indicates that measurements of R/s.<j> at the 
LHC and at the Tevatron will have strong impact on the future tuning of event generators. 
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Figure 11. Comparison of Ra<p predictions for A0 max = 5ir/Q of different event generators with 
the NLO pQCD predictions as a function of Ht in different y* regions for the LHC and Tevatron 
scenarios. 



5 Summary 

We have introduced a new quantity i?A</> Ior studies of the rapidity and transverse momen- 
tum dependence of dijet azimuthal decorrelations at hadron colliders. Since is defined 
as a ratio of cross sections, experimental and theoretical uncertainties, which are correlated 
between the numerator and the denominator, will cancel to a large extent. We suggest to 
measure the rapidity and transverse momentum dependence of i?A<f> using the longitudinally 
boost invariant variables y* and Ht- In pQCD, Ra$ is given by a ratio of three-jet and 
dijet cross sections, and, at leading order, proportional to a s . We have proposed scenarios 
for measuring R^ at the LHC and at the Tevatron, for which we have computed the NLO 
pQCD predictions and the size of non-perturbative corrections. The non-perturbative cor- 
rections are below 3% with a model dependence of typically less than 2%. The NLO pQCD 
predictions have PDF uncertainties of less than 1%, and a scale dependence of 4-6% (for 
max — 7ir/8 and 5tt/6). These properties make R/\^,, and its Ht dependence an ideal 
quantity for determinations of a s and for studying its running up to the energy frontier. 

While these studies have focused on the Ht dependence of R^, we have also investi- 
gated the theoretical predictions for the rapidity dependence. We have shown that, at fixed 
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Ht, NLO pQCD and the event generators PYTHIA and HERWIG all predict an increase of 
Ra4> with y*, however, the rate of the increase differs strongly between the different models. 
Measurements of Raa a t the LHC an the Tevatron will be able to distinguish between the 
different predictions, test NLO pQCD and play an important role in future tuning of Monte 
Carlo event generators. 
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